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Abstract 

We give an overview of recent advances in analysis of equations of 
electrodynamics with the aid of biquaternionic technique. We discuss 
both models with constant and variable coefficients, integral repre- 
sentations of solutions, a numerical method based on biquaternionic 
fundamental solutions for solving standard electromagnetic scatter- 
ing problems, relations between different operators of mathematical 
physics including the Schrodinger, the Maxwell system, the conduc- 
tivity equation and others leading to a deeper understanding of physics 
and mathematical properties of the equations. 

1 Introduction 

Application of the algebra of biquaternions to equations of electromagnetism 
has been subject of an important number of research articles and books (see, 
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e.g., i, i, m, m, m, m, m, m, m, m, m, m and many 

others). The aim of this work is to present some recent results in the field 
concerning the usage of algebraic advantages of biquaternions for analytic and 
numerical solution of Maxwell's system for chiral media as well as for inhomo- 
geneous media. Compared to a considerable number of publications dealing 
with biquaternionic reformulations of Maxwell's equations for a vacuum or 
for a homogeneous isotropic medium, application of biquaternions to elec- 
tromagnetic models corresponding to more complicated media (a much more 
challenging object for studying) was discussed in relatively few sources ([T3]. 
[21], [23], [28], [3lj, [10], HI]). Meanwhile the possibility of representation of 
Maxwell's system for a vacuum in the form of a single biquaternionic equation 
is known since 1919 [35]) only recently it became clear how this result can 
be generalized for inhomogeneous media [27], [28] and for chiral media [T3] . 
An appropriate quaternionic or biquaternionic reformulation of a first order 
system of mathematical physics opens the way for applying different methods 
which in many aspects preserve the algebraic power of complex analysis. For 
example, it is not easy to arrive at the Cauchy integral formula for holomor- 
phic functions using two-component vector formalism or even more difficult 
task to develop using this formalism a holomorphic function into a Taylor 
series. No mathematician would consider such way of presenting complex 
function theory helpful or appropriate. However, this is precisely what is 
happening in the study of three or four- dimensional models of mathematical 
physics. Compare, e.g., the Stratton-Chu integrals written in their standard 
form (see, e.g., [lO]) with their biquaternionic representation [26], [28j, ^32] 
which is in fact a convolution of a biquaternionic fundamental solution of 
the Maxwell operator with the electromagnetic field and it is quite evident 
that the latter is natural and elucidating. The meaning of the Stratton-Chu 
formulas as a Cauchy integral formula for the electromagnetic field becomes 
transparent and no doubt in their biquaternionic form the Stratton-Chu for- 
mulas could be included in a moderately advanced course of electromagnetic 
theory which is not the usual case up to now in spite of their central role in 
electrical engineering applications. 

The results presented in this work are "essentially quaternionic" in the 
sense that it is not clear how they could be obtained with other techniques. In 
the first part (sections 2-5) we explain a numerical method for solving electro- 
magnetic scattering problems with an unusual for three-dimensional models 
precision by the aid of biquaternionic fundamental solutions which to the dif- 
ference of the usually utilized matrix fundamental solutions for the Maxwell 
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equations (see, e.g., [T] and [12]) enjoy some advantageous properties. First 
of all they are not matrices but vectors (of four components). Second, they 
have a clear physical meaning of fields generated by point sources. Third, 
their singularity is lower than that of fundamental solutions based on a ma- 
trix approach. The main idea of this work is to explain how our approach 
works and how it can be used. We only formulate some necessary results like 
those about the completeness of our systems of quaternionic fundamental 
solutions in appropriate functional spaces referring the interested reader to 
some previous publications, in particular [23j where the corresponding proofs 
can be found. 

In sections 6-8 we consider the time-dependent Maxwell system for chiral 
media, rewrite it in a biquaternionic form as a single equation and then 
construct explicitly a corresponding Green function. Section 9 is dedicated 
to the time-dependent Maxwell equations for inhomogeneous media. We 
show that these equations can also be written as a single biquaternionic 
equation. In a static case the corresponding quaternionic operator factorizes 
the stationary Schrodinger operator. We study relationship between solutions 
of these important physical equations. 

2 Biquaternionic fundamental solutions 

Let E[(C) denote the set of complex quaternions (= biquaternions) . Each 
element a of EI(C) is represented in the form a = ^^=q akik where {ak} C C, 
io is the unit and {ik\ k = 1, 2, 3} are the quaternionic imaginary units: 

^0 = ^0 = -il; ioh = ikio = ik, k = 1, 2, 3; 

iii2 = -i2H = k; = -ki2 = ii] isH = -iik = ^2- 

We denote the imaginary unit in C by z as usual. By definition i commutes 
with ik, k = 0,3. 

We will use the vector representation of complex quaternions, every a G 
EI(C) is represented as follows a = + ~a , where is the scalar part of 
a: Sc(a) = Oq, and li is the vector part of a: Vec(a) = li = Yl\=i^kik- 
Complex quaternions of the form a = ~a are called purely vectorial and can 
be identified with vectors from C'^. The operator of quaternionic conjugation 
we denote by Ch'- a = = aQ — ~a . 
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Let us introduce the operator D = X]fc=i^fc^fc> where dk = gf^, whose 
action on quaternion valued functions can be represented in a vector form as 
follows 

Df = - diY f + grad /o + rot / . 

That is, Sc{Df) = - div 7 and Vec{Df) = grad/o + rot 7- 

Denote = D+a, where a is a complex constant. We have the following 
factorization of the Helmholtz operator |15j : 

A + a^ = -D^D_^ = -D_^D^. (1) 

Using the fundamental solution of the Helmholtz operator 

„ia\x\ 



(we suppose that Ima > 0), the fundamental solutions /Cq, and /C_a for the 
operators Da and D_a can be obtained from ([1]) in the following way 



]Ca = -{D-a)9a and /C_„ = -(D + (2) 

We have 

where 5 is the Dirac delta function. 

From ([2]) we obtain the explicit form of /Cq, and /C_q,: 

JC±a{x) = {±a + — ia^)6a{x). (3) 

\x\ m 

Here x = X]fc=i^A:^fc- Note that /C^ and /C_a are full biquaternions with 
Sc{]C±a{x)) = ±a9a{x) and Vec(/C-|-a(x)) = — grad^Q,(x) = {-^—ia-^^)9a{x). 

More information on the algebra of biquaternions and related calculus 
can be found in 1281. 



3 Biquaternionic reformulation of Maxwell's 
equations in chiral media 

The operators and D^a are closely related to the Maxwell equations. 
Consider the Maxwell system for a homogeneous chiral medium (see, e.g., 

mm) 
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and 



rot 'e (x) = -ia (ll (x) + piotH (a;)) (4) 
rot H (x) = ia (x) + p lot^ (x)^ , (5) 



where a = Uy/eJI. Some examples of numerical values of /3 for physical media 
can be found, e.g., in [36]. We notice only that when /3 = we obtain the 
Maxwell system for a homogeneous, isotropic achiral medium with the wave 
number ct- _^ _^ 

The vectors E and if in (jlj) and ([5]) are complex. Consider the following 
purely vectorial biquaternionic functions 

Tp = ^ + iH and ^ = E - iH . 

It is easy to verify (see [21], [31] or [23|) that Ip and ijj satisfy the following 
equations 



and 
where 



(D + ai)^ = 
(D-as) = 0, 



a a 

«! = — -rr, a2 



;i + a/3)' (l-a/5)' 



Remark 1 If P = then ai = a2 = a and we arrive at the quaternionic 
form of the Maxwell equations in the achiral case (see f3E, Sect. 9], JEE/) but 
in general ai and are different and physically characterize the propagation 
of electromagnetic waves of opposing circular polarizations. 

Obviously the vectors E and H are easily recovered from and ip : 
E = ^{^ + 7) and H = ^^{^-^). 
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4 Completeness of a system of biquaternionic 
fundamental solutions 



Let r be a sufficiently smooth closed surface in M.^. Here we use the term 
sufficiently smooth for surfaces whose smoothness allows us to introduce the 
corresponding Sobolev space if*(r) for a given s G M. 

The interior domain enclosed by F we denote by fi+ and the exterior by 

Let ~e and h be two complex vectors defined on T. 

Definition 2 We say that ~e and h are extendable into Q'^ if there exist 
such pair of vectors E and H defined on that equations and ^ are 
satisfied in fi"*" and on T we have E |r= ~e and H |r= h . 

Definition 3 The vectors ~e and h are extendable into Q~ if there exist 
such pair of vectors E and H defined on Q~ that equations and ^ are 
satisfied in Q~ , the Silver- Miiller condition 



E 



— X H 

\x\ 



oi^) (6) 
\x\ 



is fulfilled at infinity uniformly for all directions and on T: E \r= ~e and 
H \r= 'h. 

With the aid of quaternionic analysis techniques these two introduced 
classes of vector functions can be completely described. For achiral media it 
was done in [25] (see also [32i Sect. 11]) and for chiral media in [21]. Here 
we recall these results without proof. 

We will need the following operators 

Safix) = -2 j }Ca{x -y)li{y)f{y)dTy, x G F, (7) 
1 1 

Pa = -{I + Sa) and Qa = -{I-Sa) 

which are bounded in iJ^(F) for all real s. The function / in ([7j) is a biquater- 
nion valued function, it is the quaternionic representation of the outward 
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with respect to f2+ unitary normal on F: 7? = Ylt=i '"^k'ik and all the prod- 
ucts in the integrand in ([7]) are quaternionic products. From the numerous 
interesting properties of the operators Pq,, Qa and Sa (see [32]) we will need 
only the following fact 

Theorem 4 Let complex vectors ~e and h belong to H^{T), s > 0. Then 

1. in order for ~e and h to he extendable into Vf^ the following condition 
is necessary and sufficient 

{'e + it) G imP„,(iP(F)) and - it) G imP_^,(iP(F)) 

(8) 

or which is the same 

~e+ih =Sai{~e+ih) and ~e—ih = S'-ajC?— i /i ) on F. 

2. in order for ~e and h to be extendable into Q~ the following condition 
is necessary and sufficient 

Ce +ili) eimQ^^{H'{T)) and - iTT) G imQ_„,(iP(F)) 

(9) 

or which is the same 

~e+ih = —Sai{~e+ih) and ~e—ih = — S'-ajC?— ^ ^ ) on F. 

Now we will show how two systems of quaternionic fundamental solutions 
suitable for the approximation of the vector functions extendable into or 
Q~ can be constructed. 

By F~ we denote a closed surface enclosed in f2+ and being a boundary 
of a bounded domain V, and by F+ we denote a closed surface enclosing 
r2+ as shown in the figure. By {yn}^=i we denote a set of points densely 
distributed on F~, and by {y^}^]^ a set of points densely distributed on F"*". 
For each of these two sets we construct a corresponding pair of systems of 
quaternionic fundamental solutions. The pair of systems 

(10) 
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corresponds to {yn}^=i and the pair of systems 

{^a^,nix) = ^e^ii^ - yn)}'^=i and { ^Cl^^,™ (a^) = /C„„2 (x - ) }^^^ 

(11) 

corresponds to {?/„ 

The following theorems show us the possibility to apply the fundamental 
solutions ( ITOl) for the numerical solution of interior boundary value problems 
for the Maxwell equations (jl]), ([5]), and fundamental solutions ( ITTj) for the 
solution of exterior problems. 

Theorem 5 l23f Let two complex vectors ~e and h belong to H'^{T), s > 1, 
be extendable into and both a\ and be not eigenvalues of the Dirich- 
let problem in Then ~e and h can be approximated with an arbitrary 
precision (in the norm of H^~^{r)) by right linear combinations of the form 




and 



.7 = 1 



N 
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where aj and bj are constant complex quaternions. 

Theorem 6 12^1 Let two complex vectors ~e and h belong to H'^iT), s > 1, 
be extendable into and let both af and q;| be not eigenvalues of the Dirich- 
let problem in V . Then ~e and h can be approximated with an arbitrary 
precision (in the norm of H^~^{T)) by right linear combinations of the form 



^^ = ^(E^^..s + E^:..,A) (12) 



and 



"^^ = - X^^-«2,i^j)' (13) 

i=i j=i 

where aj and bj are constant complex quaternions. 

Remark 7 The right linear combinations in Theorem\^ and Theorem\^ are 
in general full quaternions. In order to ensure that they will be purely vecto- 
rial additionally to a usual boundary condition for the electromagnetic field 
we have to add the requirement that their scalar parts be equal to zero. We 
show how this can he easily achieved on some examples of numerical realiza- 
tion considered in the next section. 



5 Numerical realization 

Consider the exterior boundary value problem for the Maxwell equations 
corresponding to the model of electromagnetic scattering by a perfectly con- 
ducting body with a boundary F. Find two vectors E and H satisfying 
(jll) and (l5l) in the condition ([6]) at infinity and the following boundary 
condition 

= 7{x), X e F, (14) 

where / is a given tangential field. 

We look for the solutions in the form (fT2l) and ( JT3l) . applying the col- 
location method in order to find the coefficients aj and hj. Substitution of 
the vector part of f|T2l) in ([T^ gives us two linearly independent equations in 



E (x) X rz^(x) 
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every collocation point. In each collocation point x G F we must require also 
that 



N 



N 



(15) 



and 



-"2 J 



= 0, 



(16) 



AT N 

Sc(J]/C~^j(a;)aj - J]^- 
i=i i=i 
which gives us other two linearly independent equations. Taking into account 
that in f|T2|) and f|T3|) we have 8A^ unknown complex quantities we need 2N 
collocation points. After having solved the corresponding system of linear 
algebraic equations we obtain the coefficients aj and bj and consequently the 
approximate solution of the problem. A good approximation of the boundary 
condition (fT^ guarantees a good approximation of the electromagnetic field 
in the domain Q~ due to the following estimate (see p!2] p. 126]) 



'e 


_ + 


H 


< C 


li X E 

















where |Hloo q- stands for the supremum norm in any closed subset VL~ of Or 
and C is a positive constant depending on F and VL~ . 

The method was tested |20j, [22j, [23j using different exact solutions. For 
example, let /5 = and consequently a = ai = a2- The vectors 



and 



E'^ix) 



H'^ix) = -— rot 
ia 



rot c 9n (x) 



{X 



X G 



\{0} 



where ^ G MMs constant, represent the electromagnetic field of a magnetic 
dipole situated at the origin [ini Sect. 4.2]. They satisfy (jll) and ([5]) (for 
= 0) as well as the Silver-Miiller conditions at infinity. 
Let F be an ellipsoid described by the equalities. 



Xi = a cos 7] sin u, X2 = h sin rj sin z/, x^ = c cos 



(17) 



where < ?7 < 27r, < u </k. Then E"^ and H"^ give us the solution of the 
following boundary value problem 



rot E(x) = —iaH(x) 



X efi' 
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Tot H {x) = ia E (x) , xeQ , 



E{x) X 



n [x] 



fix), 



X eT 



where 



fix) 



C3d20aix) - 02836, 

cid^daix) - C^diO, 

C2di6aix) - 01826, 



X 



a \X J 



X 'n(x) 



We give the numerical results for a = 5, 6 = 3 and c = 2 in f[T7l) . As the 
auxiliary surface F" containing points y~ we have chosen an ellipsoid interior 
with respect to F with a, b and c multiplied by 0.15. In the following table 
we present the results for a = 1 + 0.3i and for different values of N. The 
corresponding errors represent the absolute maximum difference between the 
exact and the approximate solutions at the points on the ellipsoid exterior 
with respect to F with a, b and c multiplied by 5. 





Error for E 


Error for H 


10 


0.441E-03 


0.332E-03 


15 


0.693E-05 


0.713E-05 


20 


0.162E-05 


0.186E-05 


25 


0.245E-06 


0.248E-06 


30 


0.113E-06 


0.171E-06 


35 


0.522E-07 


0.409E-07 



A quite fast convergence of the method can be appreciated (all numerical 
results were obtained on a PC Pentium 4). 

Let us notice that the approximation by linear combinations of quater- 
nionic fundamental solutions can be applied to other classes of boundary 
value problems for the Maxwell system like for example the impedance prob- 
lem with the boundary condition 



E (x) X r?(x) 



H(x) X 'n(x) X 'n(x) = f (x) 



X e F. 



This implies some obvious changes in the matrix of coefficients of the system 
of linear algebraic equations corresponding to collocation points. 
More results and analysis of numerical experiments were given in 
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6 Time-dependent Maxwell's equations for chi- 
ral media 

Consider time-dependent Maxwell's equations 

rot E (t, x) = -dtB{t, x), (18) 

rot H{t,x) = dtD{t,x) + y{t,x), (19) 

div E{t, x) = div H{t, x) = (20) 

with the Drude-Born-Fedorov constitutive relations corresponding to the chi- 
ral media (see, e.g., [2], [33], |36j): 

5 (t, x) = fi(Hit, x) + plot H{t, x)), (21) 

D{t,x) = e{E{t,x) + f3 rot E{t,x)), (22) 

where /3 is the chirality measure of the medium. f3, e, fi are real scalars 
assumed to be constants. Note that the charge density p and the current 
density j are related by the continuity equation dtp + div j = 0. 

Incorporating the constitutive relations fl2T|) . fl22|) into the system f fTSj) - 
fl2U]) we arrive at the time-dependent Maxwell system for a homogeneous 
chiral medium 

rot H{t, x) = e{dt'E{t, x) + f3dt rot E (t, x)) + 7(t, x), (23) 
rot 'E{t, x) = -p{dtHit, x) + /3dt rot H{t, x)), (24) 

diYE{t,x) = (^^^, dwH{t,x) = 0. (25) 



Application of rot to (|23|) and fl24l) allows us to separate the equations for 
E and H and to obtain in this way the wave equations for a chiral medium 

rot rot E +epd^ E +2l3epd^ rot E +j3'^epd^ rot rot E = —pdt j —(3pdt rot j , 

(26) 

rot rot H + epd^H + 2l3epd'f rot ll + 13'^epdl rot rot ^ = rot 7- (27) 
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It should be noted that when /3 = 0, and (ITTI) reduce to the wave 
equations for non-chiral media but in general to the difference of the usual 
non-chiral wave equations their chiral generalizations represent equations of 
fourth order. 



7 Field equations in a biquaternionic form 

In this section following [13] we rewrite the field equations from Section [6] in 
a biquaternionic form. 

Let us introduce the following biquaternionic operator 

M = p^dtD + - iD (28) 

and consider the purely vectorial biquaternionic function 

V{t,x) = E{t,x) -i^H{t,x). (29) 

Proposition 8 [13] The equation 

MVit, x) = -^7 it, x) - /3^d,pit, x) + '-P^ (30) 

is equivalent to the Maxwell system the vectors E and H are so- 

lutions of ^2^) - ^2B[l if and only if the purely vectorial biquaternionic function 
V defined by /[^) is a solution of / f5^) . 

Proof. The scalar and the vector parts of fl5Ul) have the form 
-l3,/qJdtdiYE + ^divH+i{divE+f3ijdtdwH) = -I3^dtp+j, (31) 

(3./eJIdtTot^ +y^dt'E TotH -i{Tot^ +l3pdtTotH +pdtH) = -y^T- 

(32) 

The real part of fl32l) coincides with fl23l) and the imaginary part coincides 
with . Applying divergence to the equation and using the continuity 
equation gives us 

dt div H = and dt div E = -dtp. 
13 



Taking into account these two equalities we obtain from (IHTj) that the vectors 
E and H satisfy equations (^^. ■ 

It should be noted that for (3 = from (1251) we obtain the biquaternionic 
Maxwell operator for a homogeneous achiral medium for which the following 
equality is valid 

e^d^ — Ax = {^/elldt + iD){-^/elldt — iD). 

In the case under consideration {j3 ^ 0) we obtain a similar result. Let us 
denote by M* the complex conjugate operator of M: 

M* = P^/eJIdtD + + iD. 

For simplicity we consider now a sourceless situation. In this case the equa- 
tions ( l26l) and ( 1271) are homogeneous and can be represented as follows 

MM*U{t,x) = 0, 

where U stands for E or for H . 



8 Green function for the operator M 

Here we present a procedure from [13] which gives us a Green function for 
the operator M. Consider the equation 

iP^dtD + - iD)f{t, x) = 6{t, x). 

Applying the Fourier transform JF with respect to the time-variable t we 
obtain 

[(3^/ejliujD + ^Jeiiiijj — iD)E{u!, x) = S{x), 

where F{uj,x) = J-'{f{t,x)} = j^^f{t^x)e~^'^^dt. The last equation can be 
rewritten as follows 

{D + a){(3^/e]2uj — l)iF{uj, x) = 6{x), 

where a = j^^zi- The fundamental solution of Da is given by ([3]), so we 
have 

{f3y/ejlu! — l)iF{uj, x) = }Ca{x) = (a H ^ ~~ '^'^'r~^)^a{x), 
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from where 



F{uj,x) 



_{(3y/eJIuj - ly \ \x\J \xf (3y/eJIuj -1 



■I I 

llxl , , • 

47r Ixl 



We write it in a more convenient form 

1 



F{uj,x) 



where a 
A(x) 



1 



i 



c (x) 



\x\ 



E(x) 



UJ — a 

e /3 
47r|x| ' 



1-^ , 5x 



1 / ix\ X 



In order to obtain the fundamental solution f{t,x) we should apply the in- 
verse Fourier transform to F{uj,x). Among different regularizations of the 
resulting integral we should choose the one leading to a fundamental solution 
satisfying the causality principle, that is vanishing for t < 0. Such an election 
is done by introducing a small parameter y > in the following way 



fit,x) = \imJ^-'{Fiz,x)} 



(33) 



where z = uj — iy. This regularization is in agreement with the condition 
Ima > 0. We have 



:f-' {F{z,x)} = ^ r (-^ 



J A ix) + 



1 



ic{x) 



LO — a„ 



-B{x) E{x)e-^ye"^'duj 



where ay = a + iy. Expression (IMI) includes two integrals of the form 



(34) 



1 Qi^-ay ^iujt 



27r (tj - ay) 



i^duj, k = 1,2 



where c = c(x) . We have 



^ oo 

-y 



UC 



,J roo 



27T '■ ' \ ?! 

j=0 \ •> 



[yj - ay) 



j+k / ' 



(35) 
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where the change of order of integration and summation is possible because 
the two necessary for this conditions are fulfilled: the series is uniformly con- 
vergent on each segment and the integrals of partial sums converge uniformly 
with respect to j. Denote 

h,j{t) - / 

J^oo [uj — ayY 

For k = 1 and j = we obtain (see, e.g., [21 Sect. 8.7]) 

Ji,o(t) = 2TTiH{t)e'"'y 

where H is the Heaviside function. For all other cases, that is for k = 1 
and j = 1, oo and for k = 2 and j = 0, oo we have that j + k > 2 and the 
integrand in fl35l) has a pole at the point ay of order j + k. Using a result 
from the residue theory pTt Sect. 4.3] we obtain 

Ikj{t) = 2Txi Res --^ for t > and j + k > 2. 

Consider 

Res -3 = - ^ lim , ^ e-'= ^ ' , „ fort>0 

-y {oJ-ayY^" (j + A;-l)!--a,aa;^+fc-i (j + A; - 1)! 

and j + k > 2. 

For t < we have that hjit) is equal to the sum of residues with re- 
spect to singularities in the lower half-plane y < which is zero because the 
integrand is analytic there. Thus we obtain 

(jf]j+k~i 

4,(t)=2vrz/f(t)^AJ___e^«.*. 
Substitution of this result into (1351) gives us 



I,=,^H(t)e""'±L^ and I, = -H{t)e<'.H±-S^^. 

Now using the series representations of the Bessel functions Jo and Ji (see 
e.g. [m Chapter 5]) we obtain 

h = iH{t)e'''y'Jo (2Va^ and h = -H{t)^-e"'y* Ji {2^fd) . 
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Substituting these expressions in and then in (1551) we arrive at the 
following expression for /: 



f{t,x) = H{t)e'"''E (x) (^-A (x) (^2 Vet) + iB (x) Jo (2Vctj^ . 

Finally we rewrite the obtained fundamental solution of the operator M in 
explicit form: 



f{t,x 




Let us notice that / fulfills the causality principle requirement which 
guarantees that its convolution with the function from the right-hand side of 
(1301) gives us the unique physically meaningful solution of the inhomogeneous 
Maxwell system fl23l) - fl25l) in a whole space. 



9 Inhomogeneous media 

Consider Maxwell's equations in a nonchiral inhomogeneous medium. Thus 
we assume that e and n are functions of coordinates: 

e = e{x) and /i = /i(x). 

Then the Maxwell system has the following form 

rotH = e(9tE+j, (36) 

rot E = -MH, (37) 

div(£E) = p, (38) 

div(//H)=0. (39) 

In this section following the procedure exposed in we show that this 
system of equations can be written in the form of a single biquaternionic 
equation. 
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Equations fl38p and fl39p can be written as follows 

, grade _ p 
divE+ < ,E >= 

e e 



and 



divH+<^,H>=0. 



Combining these equations with (15^ and (1371) we obtain the Maxwell system 
in the form 

/)E=<^,E>-MH-^ (40) 

and 

Z^H =< ^^,H > +edt^ + l (41) 
/i 

Let us make a simple observation: the scalar product of two vectors 'j) and 
~q can be written as follows 

< ^, ^ >= -^(^M + M^)^. 

Using this fact, from (HU]) and fHT]) we obtain the pair of equations 

(Z,+ ll£|£)E = -iAf=f^E-MH-f (42) 

and 

{D + = -Im^H + .5,E + j. (43) 

Note that 

1 grad e grad i/e 



2 £ 

Then (1421) can be rewritten in the following form 



^D(v^-E) +E-7' = -/i^tH- ^, (44) 
v ^ £ 



where 

gradA/e 



e :-- 
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Analogously, takes the form 

^D{^ ■ H) + H ■ 7? = eStE + j, (45) 

where 

grad ^ 

Introducing the notations 

'S := y/eE, n := y/JIH, 
multiplying (l44l) by ^/e and (145!) by ^/JI we arrive at the equations 

{D + M^)^ = --dtH - (46) 
c 

and 

{D + M^)7^ = -^dt~S + v^TZj, (47) 

where as before c = 1/ y^e/I is the speed of propagation of electromagnetic 
waves in the medium. 

Equations fl46p and fH7|) can be rewritten in an even more elegant form. 
Consider the function 

V ■.= ~s +tn. 

Let us apply to it the biquaternionic operator 

-dt + iD. 
c 

We obtain 

{-dt + iD)\ = -dt~S - DH + i{-dtn + ot). 

c c c 

Applying (jU]) and fjl^ to the real and imaginary parts of this equation gives 



[-dt + iD)Y = -i(M^~S + iM'^n) -^i-^. (48) 
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Note that 



^ = ^(V + V*) and H = 1(V-V*). 



Hence 
Let us notice that 



^ _^ grade , _^ gradiy 

6 + fi = and e — ^ = , 

c W 

where W = ^J^j y/e is the intrinsic wave impedance of the medium. Denote 

7^:= and W := ^^^^j^ . (49) 

Then _^ 

M^^ + iM^n = -(M^V+M^V*). 

From fHHj) we obtain the Maxwell equations for an inhomogeneous medium 
in the following form 

(-dt + iDW-M'^Y-M'^V* = -(V^j + (50) 



This equation first obtained in [27], [28] is completely equivalent to the 
Maxwell system (!36 l) -( !39l) and represents Maxwell's equation for inhomo- 
geneous media in a quaternionic form. We formulate this as the following 
statement. 

Theorem 9 Let c = l/y/ejl, ~c and W he defined by (23^- Then two real- 
valued vectors E and H are solutions of the system ^SEj-^SQ^ iff the purely 
vectorial biquaternion V =y/eE + iy/JIH. is a solution of fSW . 

Note that if e and fi are constant (a homogeneous medium), equation 
fl50|) turns into a well known (at least since the work of C. Lanczos [35] ) 
biquaternionic reformulation of the Maxwell system in a vacuum which was 
rediscovered by many researchers (e.g., [E] and comments in [14j). 

Remark 10 Equation ([J^) can be considered as a generalization of the Vekua 
equation, well known in complex analysis, that describes generalized analytic 
functions fj^l - Recently in fST^ using the L. Bers approach ^ another gen- 
eralization of the Vekua equation was considered. Most likely some of the 
interesting results discussed in can be obtained for also. Their 
physical meaning would be of great interest. 
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9.1 Static case and factorization of the Schrodinger op- 
erator 



When the vectors of the electromagnetic field do not depend on time from 
(jlH]) and fH7|) we obtain two independent equations 

{D + M^)t = — ^ 

and 

Let us consider the sourceless situation, that is we are interested in the 
solutions for the operator D + M°' , where the complex quaternion "c? repre- 
sents or 7^ and has the form 



a 



grad/ 



The scalar function / is different from zero. 

Note that the study of the operator D + °' M practically reduces to that 
of D, as shown in [3H]- In the case of the operator D + M"' (which can be 
called the static Maxwell operator) the situation is quite different. 

Consider the equation 

(-A + u)g = inn (51) 

where u and g are complex valued functions, and Q is a domain in R^. We 
assume that g is twice continuously differentiable. 

Theorem 11 [29f Let f be a nonvanishing particular solution of 137]) . Then 
for any scalar twice continuously differentiable function g the following equal- 
ity holds, 

Df Df 

{D + M—){D ~ M~)g = {-A + u)g. (52) 

Remark 12 The factorization ^5M) was obtained in JB^ in a form which 
required a solution of an associated biquaternionic Riccati equation. In 
it was shown that the solution has necessarily the form Df/f with f being a 
solution of 137]) . 

Remark 13 Theorem\TT\ generalizes theorem 21 from 
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Remark 14 As g in ^5M) is a scalar function, the factorization of the Schrddinger 
operator can be written in the following form 

{D + M'T)fDif-'g) = {-A + i^)g, 



from which it is obvious that if g is a solution of l[51\) then the vector F = 
fD{f~^g) is a solution of the equation 

Df 

{D + M~)F = mVL. (53) 

The inverse result is given by the next statement where the following notation 
is used 

X y z 

A[G]{x,y,z) = jGi{^,yo,zo)d^ + J G2ix, (, Zo)dC + j G3{x,y,r])dr] + C 



yo zo 



(C is an arbitrary complex constant). 

Theorem 15 J^W Let F be a solution of in a simply connected domain 
Q. Then g = /^[/~^F] is a solution of 137]). 

Moreover, a factorization of the operator divpgrad+g is valid also. 

Theorem 16 Let Uq be a nonvanishing particular solution of the equa- 
tion 

(divpgrad+g)M = inVlcM.^ (54) 

with p, q and u being complex valued functions, p G C^(f2) and p in Q. 
Then for any scalar function ip G C'^{Q) the following equality holds 

(divpgrad+g)^ = -p^^^{D + M—){D - M~)p^'^^ (55) 

where f = p^^'^Uq. 

Thus, if M is a solution of equation fl5^ then 

F = fD{f-y/\) = fD{u,\) 

is a solution of equation fl53l) (see remark . The inverse result has the 
following form. 
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Theorem 17 J^P| / Let F be a solution of equation in a simply connected 
domain f2, where f = p^^^Uq and Uq be a nonvanishing particular solution of 
3. Then 

u = uoA[f-^¥] 

is a solution of 

Notice that due to the fact that in (l55l) 99 is scalar, we can rewrite the 
equahty in the form 

(divpgrad+g)^ = -p^''^{D + M^){D - ^Ch)p^''^'^. 

Now, consider the equation 

Df 

(D - -jCh)W = 0, (56) 

where W is an E[(C)-valued function. Equation (156|) is a direct generahzation 
of the main Vekua equation considered in [30]. Moreover, we show that it 
preserves some of its important properties. 



Theorem 18 1291 Let W = Wq + be a solution of ([5|j. Then Wq is a 
solution of the stationary Schrddinger equation 

- AWo + uWo = 0, (57) 

where v = A///. Moreover, the function u = f^^Wo is a solution of the 
equation 

div(/2gradu) = (58) 
and the vector function v = /W is a solution of the equation 

rot(/-2rotv) = 0. (59) 

Remark 19 f2^ Observe that the functions 

p-f p-^i p_^3 
-TO — /, -Ti — y, -^2 — y, ^ ~ Y 



give us a generating quartet for the equation f56]l . They are solutions of (l5i 
and obviously any M.{C) -valued function W can be represented in the form 



3 

j=0 
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where ipj are complex valued functions. It is easy to verify that the function 
W is a solution of / f55]) iff 

3 

J2 (DVj) F, = (60) 

j=0 

in a complete analogy with the two-dimensional case. Denote 

Then [6D\) can he written as follows 

D{w + w)f + D{w - lZJ)i = 
which is equivalent to the equation 

1-/' 

Dw = 

Remark 20 The results of this section remain valid in the n-dimensional sit- 
uation if instead of quaternions the Clifford algebra C/o,n (see, e.g., JS^, fWj) 
is considered. The operator D is then introduced as follows D = Yl^=i ^igfr 
where Cj are the basic basis elements of the Clifford algebra. 

Acknowledgement 21 The authors wish to express their gratitude to CONA- 
CYT for supporting this work via the research project 50424- 
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